Asperity contacts at the nanoscale: comparison of Ru and Au 
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We develop and validate an interatomic potential for ruthenium based on the embedded atom 
method framework with the Finnis/Sinclair representation. We confirm that the new potential 
yields a stable hep lattice with reasonable lattice and elastic constants and surface and stacking 
fault energies. We employ molecular dynamics simulations to bring two surfaces together; one flat 
and the other with a single asperity. We compare the process of asperity contact formation and 
breaking in Au and Ru, two materials currently in use in micro electro mechanical system switches. 
While Au is very ductile at 150 and 300 K, Ru shows considerably less plasticity at 300 and 600 
K (approximately the same homologous temperature). In Au, the asperity necks down to a single 
atom thick bridge at separation. While similar necking occurs in Ru at 600 K, it is much more 
limited than in Au. On the other hand, at 300 K, Ru breaks by a much more brittle process of 
fracture/decohesion with limited plastic deformation. 
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I. INTRODUCTION 



Contacts between surfaces are central to a wide range 
of technologies and present an interesting array of fun- 
damental issues. Contact behavior may involve phenom- 
ena as diverse as adhesion, friction, bonding, wear and 
fracture. Most contact problems are intrinsically mul- 
tiscale; involving atomic bonding, nano-scale asperities, 
defect nucleation, elastic and plastic deformation, sur- 
face roughness, far- field loading, etc. The atomistic or 
nano-scale regime is particularly relevant since this is the 
length scale at which the first contact between rough sur- 
faces occurs. Furthermore, recent technological applica- 
tions like Micro Electro Mechanical Systems^ (MEMS) 
use mechanical contacts that are only a few micrometers 
large, such that nanoscale contact physics plays an even 
greater role. 

Because of its large electrical conductivity, gold has 
become a popular material for electrical contacts. Nev- 
ertheless, its application to micro-contact technology like 
MEMS switches has been hindered by a lack of reliabil- 
ity. Extensive surface damage leads to failure of some 
gold contact devices after only several million open/close 
cycles.^ The resulting surface damage has been studied 
experimentally with atomic force microscopes in order 
to understand the effect of adhesion, thermal dissipation 
and contamination-^^J^ Alloying Au with other metals 
results in increased hardness, b ut also in increased resis- 
tivity.'^ At omic- level simul at jpn j '^^ ^'^' ^^i^ and experimen- 
tal observations^^ESI j^ayg shown that the separation of 
gold contacts is accompanied by substantial plastic defor- 
mation, leading to ductile material separation with con- 
siderable material transfer from one side of the contact 
to the other. This results in significant contact surface 
morphology evolution in each cycle. Other pure metals, 
like ruthenium, have recently been used to build metal 
contacts in MEMS switches. Ruthenium contacts have 
proven to be more reliable, routinely surviving millions 
of cycles without significant degradation of the contacts.^ 
While Au contacts have been the focus of several atom- 



istic simulations, to our knowledge, there have been no 
atomistic simulation studies of contact between Ru sur- 
faces. Therefore, we do not know whether the fundamen- 
tal contact evolution mechanisms observed in Au occur 
for other metals generally, and to the important case of 
Ru, in particular. Information of this type can help guide 
the development of future generations of MEMS contact 
switches. 

A comparison of Au and Ru contact behavior is also in- 
teresting because the crystal structures of Ru (hexagonal 
closed packed) and Au (face centered cubic) differ: this 
can lead to fundamental differences in crystal plasticity. 
Molecular dynamics (MD) simulation is a powerful tool 
for studying adhesion, defect formation and deformation 
on the nano-scale level (that of asperities in MEMS con- 
tacts). The first step in performing large scale molecular 
dynamics simulation of Ru is to establish an appropri- 
ate interatomic potential. We have investigated several 
potentials for Ru available from the literature, but en- 
countered problems with each of sufficient magnitude to 
make them inappropriate for contact simulations. In this 
article, we develop a new embedded atom method^^*^ 
(EAM) type potential for Ru. We implement this po- 
tential in our molecular dynamics simulations of asperity 
contact and separation. The main focus of this paper is 
a comparison of how asperity contacts form and separate 
in gold and in ruthenium. 



II. SIMULATION METHOD 

In order to understand the formation and separation 
of nanoscale asperity con tacts, we employ the molecular 
dynamics (MD) technique^^^ to simulate the accompa- 
nying atom-scale dynamics. In MD, the motion of the 
atoms are simulated by solving Newton's equations of 
motion for each atom in the simulation cell. The molec- 
ular dynamics were carried out at a constant pressure and 
temperature. This was done using the Nose-Hoover ther- 
mostat!^ and barostatP^with a time step of ^t=0.0025 ps. 
In the contact experiments, we simulate two substrates. 



of roughly 15 close-packed planes, facing one another. 
These correspond to the (0001) planes of the Ru hexag- 
onal close packed (hep) lattice and to (111) planes of the 
Au face centered cubic (fee) lattice. On one of the sub- 
strates we place an asperity, as shown in Fig. [l] The top 
two layers in the upper substrate and the bottom two in 
the lower substrate were kept rigid in order to prevent 
the two substrates from rotating and to fix the relative 
displacement of the two substrates in the z-direction. Pe- 
riodic boundary conditions are enforced in the x and y- 
directions. The simulation cell contains 35,000 atoms. 
In the present simulations, an anisotropic isobaric ensem- 
ble is used at zero lateral pressure to maintain constant 
(zero) stress in the x and y directions. These boundary 
conditions were chosen to model a periodic surface while 
allowing the system to expand or contract in the x and 
^/-directions, as it would if the substrates were of finite 
extent in these directions. 

Our first step in preparing the system is to construct 
a cube-shaped asperity with 8 close packed layers on the 
upper surface of the bottom substrate. Second, we anneal 
the asperity by increasing the temperature from 3 K to 
1700 K in 100 ps. Next, we quench the system back to 3 
K within 100 ps, and finally we increase the temperature 
to 600 K or 300K in 100 ps, for Ru and Au, respectively. 
The system is then quenched to the target temperature 
that may differ from the above values. This approach 
allows the asperity shape to evolve towards one that is 
relaxed away from its initial cubical shape. The resultant 
Ru asperity at 300 K is shown in Fig. [l] 
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FIG. 1: (Color Online) Atomic configuration of the simula- 
tion cell after annealing of the asperity to a final temperature 
T=300 K and prior to the contact simulation experiment. The 
top and bottom substrates consist of the same material but 
different colors are used to emphasize material transfer after 
contact separation. 

After the asperity annealing is completed, we displace 
the upper substrate toward the lower substrate at a con- 
stant velocity of 0.07 A/ps (7 m/s), while holding the 
lower substrate fixed (i.e., we hold the two atomic lay- 
ers at the bottom of the lower substrate fixed). When 
the distance between the top of the upper substrate and 



the bottom of the lower substrate is 1.37 nm and 1.23 
nm, for Ru and Au, respectively, the sign of the veloc- 
ity of the upper substrate is reversed. The simulation 
continues until the upper and lower substrates are com- 
pletely separated. The substrate velocity used here corre- 
sponds to that in a MEMS switch operating at :^ 1 GHz. 
The 2:-component of the force on the upper substrate is 
calculated during the simulation in order to determine 
the force-displacement relation for the contact simula- 
tion system. We investigate the formation of defects us- 
ing the order parameter recently developed by Ackland 
and Jones'^. 



III. INTERATOMIC POTENTIALS 



The gold interatomic potential employed in the present 
study is an embedded atom method (EAM) type poten- 
tial developed by Cai and Ye^. This potential provides 
a reasonable description of both bulk and surface phe- 
nomena an d was u sed previously in contact simulation 
exper iment s .'^'^'^tlS 

We examined several potentials for ruthenium that 
were available in the literature. For the present simu- 
lations, we were looking for a potential that reproduces 
the elastic constants and cohesive energy of Ru, and also 
gives stacking fault energy and surface energies that are 
in reasonable agreement with experiment. Stacking fault 
energies play an important role in plastic deformation 
and surface energy determines the work of adhesion - 
both key elements in contact formation and separation 
phenomena. Igarashi et al.^^ developed two different 
EAM (Finnis-Sinclair) potentials for Ru. The stacking 
fault energies were /2= 13 [meV/A^] in one version and 
/2= 30.54 [meV/A^] in the other. Both of these dif- 
fer significantly from the experimental value of /2= 53.8 
[meV/A^]. Hence, this potential is unsuitable for our 
contact simulations. While modified embedded atom 
method (MEAM) potentials'^ '^ have significantly more 
fiexibility, the MEAM potentials fail to produce a stable 
hep crystal, as noted previously by Mae et al.'^. We con- 
firmed this via MD simulations, which showed that appli- 
cation of a small strain destabilized the hexagonal close 
packed lattice of MEAM Ru. Grinberg et al. '^ developed 
an EAM potential for Ru, but did not include (or calcu- 
late) the surface and stacking fault energies. Our expe- 
rience shows that such properties will not be reproduced 
unless they (or surrogate properties) are included in the 
fit. Hu et al.^'' recently proposed another Ru potential, 
but the calculated values of the stacking fault energy /2 = 
4.25 [meV/A'l and the surface energy of the basal plane 
7= 80 [meV/A^] are far from the experimental values. 

Because none of the extant Ru potentials met the re- 
quirements for contact simulations, we develop a new 
EAM potential that reproduces the stacking fault energy 
and the surface energy of ruthenium with good accuracy. 



TABLE I: The analytical form of the Ru potential. All distances are expressed in A and energies in eV. The two numbers 
listed as "cutoff' indicate the range of the adjacent basis function: a function with cutoff x — y should be multiplied by 
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where H is the Heaviside step function. The poteni 
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Function 


Value 


Cutoff 


V(r) 


exp(11.393523549007 - 8.8769414996354r + 3.7751688938761r2 - 0.99892512021332r2) 


1.0-1.9 




+84.370697642971(2.9 - r)^ - 370.62682398816(2.9 - 


-rf 


1.9-2.9 




+655.67905839695(2.9 - rf - 549.25562909352(2.9 - 


-rf 


1.9-2.9 




+177.09076987937(2.9 - r)^ 




1.9-2.9 




-3.0363100686117(5.0 - r)^ - 3.0909218331545(5.0 - 


-rf 


1.9-5.0 




-4.4211545804424(5.0 - r)^ - 1.9205251191921(5.0 - 


-ry 


1.9-5.0 




-0.31060214601543(5.0 -r)^ 




1.9-5.0 




+0.93316319033507(6.4 - r)^ - 3.0761272816964(6.4 


-rf 


1.9-6.4 




+3.5391800334269(6.4 - r)^ - 1.7626076559210(6.4 - 


-ry 


1.9-6.4 




+0.32587864688719(6.4 - r)^ 




1.9-6.4 


$(r) 


1.180972686(4.8 -r)^ 




0-4.8 




-0.7740233148(4.8 - r)^ + 0.2696095863(4.8 - rf 




0-4.8 




-0.0471950464(4.8 - rf + 0.0032605667(4.8 - rf 




0-4.8 


F(p) 


0.88124070590302p°-^ 




O-oo 




+6.5448285611211 x 10-^(p - 65)^ 




65-00 




-2.5046137745775 x 10"^(p - 75)^ 




75-00 




+2.3771143627356 x 10"^(p - 90)^ 




90-00 




-4.4929936875438 x 10"^ (p - 95)^ 




95-00 




+2.3927653989776 x 10"^(p - 100)^ 
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The EAM potential between the atoms is given b; 



N-l N 



N 



^=E E v{n,)^^F{p,) 



(1) 



= 1 i=i+i 



where the subscripts i and j indicate each of the N atoms 
in the system, r^j is the distance between atoms i and 
j, V{rij) is a pairwise potential, F{pi) is the embed- 
ding energy function and pi = ^a ^{^ij) is the radially- 
symmetric electron density at the position of atom i. The 
function <l>(r) is another pairwise potential. Creating 
the new potential involved finding optimal functions for 
V{rij), F{pi), and ^{rij) of Eq. ([T]). The potential de- 
velopment procedure was described in details in Ref. [29l 
The analytical form of the functions for Ru is reported 
in Table [11 

Table [ll] shows a comparison between several crystal 
properties determined from the potential and values ei- 
ther determined from calculations or experiments. Note 
that lattice constants, surface energy 7 and basal stack- 
ing fault energy I2 are in ver y good agreement with 
the experimental values . ^ l ^^ * ^^ * We obtained larger dif- 
ferences between our calculated values and experimen- 
tal values of the other properties. In particular, we 
determine the melting point of this Ru potential to be 
Tjn = 1792 + 5 K, using the coexistence method. ^^ While 
this is low compared with the experimental value for Ru 
Tjn = 2607 K, this is not a problem since all simula- 
tions were performed at temperatures less than a quar- 



ter of the melting temperature. Overall the potential is 
much better suited for contact simulations than previ- 
ous potentials. Furthermore, we performed a series of 
MD simulations using this potential over a wide range of 
temperatures in order to insure that the hep structure 
is stable at finite temperatures. These simulations were 
performed on a perfect crystal containing 4800 Ru atoms 
and using periodic boundary conditions in all three di- 
rections, where the periodic lengths were free to adjust 
to insure zero stress in the system. The hep structure 
was found to be stable from low temperature up to the 
melting point. 



IV. RUTHENIUM CONTACTS 

Figures [2] shows the force-displacement curves for the 
Ru contact simulations at temperature T=300 K and 
T=600 K. The two force- displacement curves have many 
features in common. The zero of the displacement cor- 
responds to the first contact between the top substrate 
and the asperity; at large separation (initially negative 
displacement), the force is zero. As the two substrates 
approach each other, the force becomes negative because 
of the short range attractive interactions between the 
atoms on the opposing substrate surfaces (this is deter- 
mined by the interaction range of the interatomic po- 
tentials). This attraction elastically str etches t he two 
materials and there is a jump-to-contact.'^^'SSMI ^s the 



TABLE II: Values of several physical properties of ruthenium determined using the new potential developed herein compared 
with values from calculations or experiments. The properties include the lattice constants a and c (Ref. 30), elastic constants 
(Cii, Ci2, Ci3, C33, C44, Ref. 31), (0001) surface energy 7 (Ref. _32J, cohesive energy Ec (Ref. 33J, vacancy formation en- 
ergy EJ (Ref. 22), basal stacking fault energy I2 (Ref. |34|, energy difference between the hep and fee lattices AEhcp->fcc, 
energy di fference between the h ep and bcc (body centered cubic) lattices AEhcp->bcc (Ref. i35i) . and melting temperature Tm 
{ http://www.webelements.coml) . All of these properties were used in fitting the potential, except for the melting temperature. 





Literature T=0 K 


T=3 K 


T=300 K 


a [A] 


2.706 


2.705 


2.7046 


2.7067 


c[A] 


4.282 


4.288 


4.288 


4.288 


c/a [A] 


1.582 


1.585 


1.5854 


1.5841 


I2 [meV/A^] 


54.6 


53.8 






7 [meV/A^] 


190 


189 


189 


189 


Cii 


563 


552^ 


514 


476 


C12 


188 


165^ 


148 


152 


Ci3 


168 


170^ 


145" 


162" 


C33 


624 


611^ 


586 


610 


C44 


181 


198^ 


199 


197 


Ec [eV] 


-6.74 


-6.864 


-6.864 


-6.777 


EJ [eV] 


1.85 


2.218 






AEhcp->fcc [eV] 


0.125 


0.152 


0.153 


0.137 


AEhcp->bcc [eV] 


0.265 


0.216 


bcc unstable t 


cc unstable 


Tm [K] 


2607 


1792 







(a) Average of C13, C31, C23- 

(b) Calculated using a virial expression. 



displacement increases beyond this point, the material 
goes into compression. There is an approximately lin- 
ear rise in the force with displacement, punctuated by a 
series of relatively sharp drops. The linear increase be- 
tween the drops corresponds to elastic compression. The 
sharp drops correspond to defect generation or annihi- 
lation events (see below). When the sign of the veloc- 
ity changes, the system begins to unload from the com- 
pressed state. The unloading is initially characterized by 
linear elastic regions with a few jumps. Eventually, the 
force reaches a minimum (i.e., a maximum tensile stress), 
following which the tensile stress slowly decreases to zero 
over a long displacement range. In this region, the over- 
all force-displacement trend is also interrupted by sharp 
jumps. 

The two substrates separate at larger (negative) dis- 
placement than that at which the initial contact occurred 
(i.e., zero displacement). The magnitude of the displace- 
ment necessary to separate the two substrates is much 
larger in the high temperature simulation than for the low 
temperature one, suggesting that the material is much 
more ductile (larger strain to failure) at high temper- 
ature than at low temperature. This interpretation is 
easily confirmed by reference to the morphologies ob- 
served during deformation at 300 K and 600 K shown 
in Fig.[3ta)-(f) and Fig.[3];g)-(1), respectively. At T=300 
K (Fig!~^a)-(f) ) we observe a fracture-like rupture of 
the bridge of atoms connecting the two substrates. On 
the other hand, at T=600 K (Fig. p[g)-(l)) we observe 



the formation and plastic stretching of a relatively long 
neck between the two substrates. 

Before we begin the examination of the atomic config- 
uration and how it evolves during contact and contact 
separation, we first remind the reader of the hep struc- 
ture, as viewed in the projection employed in many of the 
images in the remainder of the paper. Figure |4] shows the 
hep structure viewed along the [1120] direction (i.e., the 
(1120) plane). In this view, basal (0001) and pyramidal 
(1011) planes along which dislocations can glide are visi- 
ble. We will also show views along the [0001] directions; 
i.e., the basal planes. The [0001] projection shows dislo- 
cations thanks to the atom mismatch in different planes. 

To determine the microscopic origin of the force drops 
in the force-displacement curve, we examine the structure 
for the presence, formation and annihilation of defects. 
Based on earlier contact simulation observations^ ^^ and 
the fact that dislocations carry plastic deformation, we 
focus upon dislocations. The local atomic environment 
around dislocations and the stacking faults that they can 
create differ from those in the perfect crystal. Therefore, 
we employ the order parameter developed by Ackland 
and Jones ^^ to describe the environment around individ- 
ual atoms since it is able to distinguish between different 
local configurations (including those of different crystal 
structures). The presence of local structures of differ- 
ent or unknown symmetry within the metal indicates the 
presence of a defect. In the images discussed below, we 
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FIG. 2: (Color Online) Force-displacement diagram for Ru 
contacts for (a) T=300 K and (b) T=600 K. Positive (nega- 
tive) force indicates compression (tension). The capital letter 
labels in the plots correspond to configurations discussed in 
the text. The arrows indicate the direction in which the curves 
are traversed. 



color code the atoms according to this symmetry param- 
eter, coloring atoms red in an hep environment, blue in 
an fee environment and yellow, otherwise. A stacking 
fault along the (0001) planes of the hep structure will 
appear as a plane of blue atoms in the red matrix. Dis- 
location cores and stacking faults in other planes (such 
as the pyramidal planes) will appear yellow. 

Figure |5|a) shows the atomic configuration in a four 
(1120) atomic plane thick slap which cut the system near 
its center in a direction orthogonal to the substrate. This 
configuration corresponds to the point A in the force- 
displacement curve of Fig. |2|a). Figure [5|b) shows the 
projection of the system at the same point on a (0001) 
plane parallel to the substrate. When stress is applied 
in the [0001] direction, we find that dislocations nucleate 
in the corners of the asperity /substrate contact, as seen 
in Fig. [5] (a) and (b). Careful examination of the dislo- 
cation produced shows that the active slip system here 
is {1011} (1213) indicating the presence of dislocations 



with Burgers vector of the|(1123) type. These disloca- 
tions dissociate on the pyramidal and basal plane .^ ^^*^^ l 
The small force drops in the force-displacement curve are 
associated with plastic events in the system, like the for- 
mation or annihilation of these dislocations. 

Figure |5|c) and |5|d) show the configuration corre- 
sponding to the point B in Fig. [2|a). The large force 
drop at B occurs due to the disappearance of half of 
a (0001) plane at the interface - originally as seen in 
Figs. [5|a),(c). The edge dislocation in the basal plane 
disappears at point B [Fig. |5|c),(d)], leaving behind an 
asperity with one less (0001) atomic plane. We observe 
that, despite the disappearance of this plane, defects ac- 
cumulate within the upper substrate. Careful examina- 
tion of Fig. [5|d) shows that at least four lines of yellow 
atoms (dislocations) on different {1011} planes are vis- 
ible. This is not surprising since several {1011} planes 
have the same orientation with respect to the loading 
axis. This results in an asperity that is broadening sym- 
metrically during loading. The large force drops during 
loading are associated with the disappearance of {0001} 
planes from the asperity. (Similarly, large force jumps 
during contact separation correlate with the creation of 
new {0001} atomic planes in the asperity.) 

During loading, defects continue to accumulate in the 
proximity of the asperity as shown in Figs. [6Fa),(b)). 
This state corresponds to point C in Fig. [2|a). Between 
point C and D, a force drop occurs, corresponding to the 
decrease of one atomic plane in the asperity as shown in 
Figs. |6|a),(c). This takes place through a series of slip 
events that have the effect of expelling the defects outside 
of the asperity region. (Figs. [6]^c),(d)). 

The processes shown in Figs. [5] and [6] repeat until the 
contacts begin to separate. In the early stages of the un- 
loading process, the force decreases linearly as the top 
substrate is retracted. At T=300 K, when the unload- 
ing process is started, dislocations are still present in the 
system. These defects survive the early stages of the 
unloading and at point E in Fig. [2|a), the dislocations 
are extended, as shown in Fig. [7Fa),(b). At point F in 
Fig.[2|a), the defects have almost completely disappeared 
through a series of slip events in the pyramidal planes 
(Fig. l7Fb),(c)). During this process, a new plane is cre- 
ated inside the bridge of atoms, and a large force jump is 
observed in the force-displcement curve. After this first 
event, the elastic response continues until the slope of the 
force-displacement curve changes. At this point, we ob- 
serve both crack propagation along with slip. Both cracks 
and dislocations are nucleated at the asperity/substrate 
contact (Fig. Isl). In particular, from the top view, we 
observe how plastic events occur mainly near the outer 
surface of the asperity, while in the center we observe 
fracture, as seen from the near perfection in the cross- 
sections in Fig. |8Fb),(d). These fracture events correlate 
with small jumps in the force-displacement curve. 

When the contact simulation is performed at T=600 
K, we observe the same general features as seen at the 
lower temperature, T=300 K, during the loading portion 
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FIG. 3: (Color Online) A series of configurations corresponding to decreasing displacements d in Angstroms, during the 
separation of the Ru contact surfaces at (a)-(f) 300 K and (g)-(l) 600 K. The atoms are all Ru, but are colored red or yellow 
to indicate that they were initially part of the upper or lower substrates, respectively. 
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FIG. 4: (Color Online) Schematic representation of the [1120] 
projection of the hep structure. The single and double red 
lines indicate a (0001) basal plane and a pair of 1011 pyra- 
midal planes, respectively. The letters, A and B, refer to 
stacking sequence of the basal planes. The coloring is imply 
to distinguish between these planes. 




of the cycle. The same slip events occur, but with less 
accumulation of defects in the substrate. This is most 
likely the result of the larger dislocation mobility at el- 
evated temperature. On the other hand, the unloading 
process is very different. Although the same slip system, 
{1011} (1213), is active at T=600K, the deformation at 
T=600K is not accompanied by fracture. Dislocations 
are nucleated at the corner of the bridge/substrate con- 
tact. The bridge elongates by repeated slips along these 
planes and the rearrangement of the atoms to form new 
atom layers (plastic accommodation of the tensile strain 
in the z-direction). An example of the slip process is 
shown in Fig. [9J 



FIG. 5: (Color Online) Atomistic configuration during the 
compression of a Ru asperity (lower substrate) against a flat 
Ru surface at T=300 K. Red atoms are in an hep environment, 
blue atoms within an fee environment and yellow atoms are 
in an environment of some other, unknown symmetry, (a) 
and (b) correspond to point A in Fig.lSJa). (a) Lateral view 
of a thin slab of four (1120) planes, while (b) the top view of 
the contacts projected on the basal (0001) plane at point A 
in Fig.[2[a). (c) Lateral view of a thin slab of 4 (1120) planes 
at point B in Fig. |2ra). (d) The top view of the contacts 
projected on the basal (0001) plane at point B in Fig. [21 a). 
Lines in (a),(c) are guides to the eye. The arrows in (bj,(d) 
indicate the position of the slab section reported in (a),(c). 



V. GOLD CONTACTS 

In order to compare the behavior of Ru to that of the 
well studied case of Au, we repeat the contact simulations 







FIG. 6: (Color Online) Same as figure ^ but for point C 
(a),(b) and point D (c),(d) in figure l2|a). 



FIG. 8: (Color Online) Same 
at point G (a),(b) and point 



as Fig. [5] but durii 
H(c),(dTinFig.[2l; 



during unloading 
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FIG. 7: (Color Online) Same as Fig.js] but during unloading, 
at point E (a),(b) and point F (c),(d) in Fig. |2ja). 



for Au at approximately the same homologous tempera- 
ture T/i = T /Tm^ where T^ is the melting temperature. 
Since, the melting temperature of Ru is T^=2607 K and 
Au melts at T^=1337 K we now analyze Au at T=150 K 
and T =300 K. To make the conditions as similar as pos- 
sible between Au and Ru we apply the stress in the {111} 
direction of Au, that is perpendicular to the hexagonal 
plane of the fee lattice which is akin to the basal {0001} 
plane of hep Ru. Since the effect of orientation was stud- 
ied in detail before^, we only briefly analyze the mor- 
phology and force-displacement curve of {lll}-oriented 
Au here. 

The force displacement curves for {lll}-oriented Au at 
T=150K and T=300K are shown in Figs. ^ The long 
tails during the separation process indicate that ductile 




FIG. 9: (Color Online) Same as Fig. [5] but for temperature 
T=600 K during unloading. Points A and B in Fig. [2lb) are 
shown in (a),(b) and (c),(d), respectively. 




FIG. 10: (Color Online) Force- displacement diagram for Au 
contacts at (a) T=300 K and (b) T=150 K. 





to accommodate strain along any axis (see Fig. [T3[a)). 
(a) (b) 



FIG. 11: (Color Online) Simulation images of Au contact 
separation at (a) T=300 K and (b) T=150 K. Red (yellow) 
spheres indicate gold atoms that at the beginning of the sim- 
ulation are in the top (bottom) substrate. 



separation (without crack propagation) occurs at both 
temperatures. This is confirmed by observation of the 

A long, thin, connective neck 



morphology in Fig. 11 



forms and elongates as the two substrates are pulled 
apart. The neck breaks when it thins to a diameter of 
1-2 atoms. This kind of separation process is responsible 
for the observed long tail in the force-displacement un- 
loading curve for Au. The amount of material transfered 
between the two surfaces is about the same at the two 
different temperatures. Since gold has an fee lattice, par- 
tial dislocations tend to form and move on {111} planes. 
Indeed, we find that with the stress oriented in the [111] 
direction, the microscopic mechanism for plas tic d efor- 
mation is the same as was previously observed^ ^^ with 
the stress in the [100] direction. 



VI. DISCUSSION 

Operationally, we consider a material to be brittle 
when separation occurs without plasticity and ductile 
when there is considerable plasticity and fracture/bond 
breaking is very limited. Figures [12] (a) and (b) show 
schematic representations of brittle and ductile separa- 
tion, respectively. Clearly, different degrees of ductility 
are possible, depending on the number of activate slip 
systems and the dislocation mobility at the test temper- 
ature. Materials that are macroscopically brittle may un- 
dergo some plasticity on the microscopic level.^^ In bulk 
metals, plasticity is normally controlled by the motion 
of dislocations. However, in nanoscale metals, plastic- 
ity is often controlled by the nucleation of dislocations, 
rather than by their motion. Gold typical exhibits duc- 
tile behavior j -^^^ l ^^ l ^^ l with necking at large strains at both 
T=150 K and T=300 K. Our simulations show that while 
Ru is ductile at T=600 K, it is more brittle at T=300 
K, separating by a combination of fracture and plastic- 
ity. The difference between the ductile/brittle behav- 
ior of Au and Ru contact has also been seen in finite 
element method calculations.^^ There are several con- 
tributions to the difference in behavior of Au and Ru. 
Since Au is fee and deforms by dislocation motion of the 
{111} (Oil) type, there are many slip systems available 






FIG. 12: (Color Online) Schematic representaton of the (a) 
brittle (several intersecting cracks) and (b) ductile separation 



scenarios. 




FIG. 13: (Color Online) (a) Two {111} slip planes in the fee 
(planes A and B) lattice, (b) Several different possible slip 
planes in the hep lattice; A is the (0001) slip plane, B is the 
(0111) slip plane, C is the (1122) slip plane, and D is the 
(0110) slip plane. 

On the other hand Ru has an hep lattice, and although 
a number o f differe nt slip systems have been observed 
(Fig. [l3];b)),l2SElESI some of these are difficult to acti- 
vate. As a result, the total number of slip systems that 
participate in deformation can be activated at low stress 
is more limited than in fee. We have also performed a 
series of simulations in which a Ru contact was loaded 
along the [1012] direction such that both the basal and 
pyramidal plane are active slip systems. Here too we ob- 
served the Ru contact to separate in a brittle way at room 
temperature and with more ductility at T=600 K. The 
surface energy of Ru is approximately 3 times larger than 
in Au and the stacking fault energy of Ru is roughly 16 
times higher than in Au.^^ Increasing the stacking fault 
energy makes slip more difficult, while increasing the sur- 
face energy (and work of adhesion) makes fracture more 
difficult. However, since the increase in stacking fault 
energy in Ru comparatively to Au is much greater then 
the increase in surface energy, the change from Au to Ru 
leads to an increased tendency for fracture in Ru than 
in Au. We speculate the increased resistance of Ru to 
plastic deformation during contact may contribute to its 
increased resistance to damage during repetitive contact 
(i.e., the morphology changes little on cycling in Ru, un- 



like Au). 



VII. CONCLUSIONS 



We have investigated the nano-contact mechanics of 
ruthenium using molecular dynamics simulations. We 
first developed a new EAM potential for Ru that re- 
produces the experimental elastic constants, surface, and 
stacking fault energies. The simulations were performed 
in a system composed of two substrates, one with a flat 
surface and the other with a parallel flat surface that has 
a single, roughly hemispherical, asperity. During load- 
ing, the top flat substrate moves toward the lower sub- 
strate at a constant velocity. At a preset displacement, 
the sign of the velocity is switched and the two surfaces 
pull apart, until the contact is broken. During the sim- 
ulation, the force in the system, the atomistic configu- 
ration, and the local order parameter are monitored. In 
Ru at T=300 K and T=600 K, we observed plastic events 
(dislocation formation, motion, and annihilation/escape) 
during the loading process. This deformation was local- 
ized primarily as slip on pyramidal planes. When the 
two substrates are pulled apart, we observed brittle con- 
tact separation at T=300 K. This was characterized by 
crack formation and bond breaking and by a short tail in 
the force-displacement curve past the maximum tensile 
force. On the other hand, at T=600 K, the separation 
was much more ductile and was characterized by slip on 
pyramidal planes and no cracks. A small bridge of atoms 
between the top and the bottom substrates forms in the 



final stages of the separation at T=600 K. The force- 
displacement curves showed considerable plastic defor- 
mation following the peak tensile force - associated with 
the necking that led to neck formation. We repeated 
nearly identical contact simulation in Au at the same 
homologous temperatures, for the sake of comparison. 
At both T=150 K and T=300 K, Au is strongly ductile 
with a long tail in the force-displacement curve and con- 
tact separation occurs with the formation of a neck. In 
Au, the plastic deformation corresponds to slip along the 
{111} planes. Our findings provide a partial explanation 
of why Ru contacts are more reliable than Au contacts. 
That is, Ru undergoes much less plastic deformation and 
morphology change than does Au under the same condi- 
tions. To be a good contact, the material must have a 
high electrical conductance. Unfortunately, the conduc- 
tance of Ru and its failure are sensitive to oxygen and 
carbon contamination.^^ In real Ru contacts, the sur- 
faces are typically covered with layers of ruthenium oxide 
that can change the adhesion and stacking fault energies. 
Work is in progress to develop a model to study the ef- 
fect of ruthenium oxide on the mechanical properties of 
nanocontacts. 
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